Transient dynamics and momentum redistribution in cold atoms 
via recoil-induced resonances 
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We use an optically dense, anisotropic magneto-optical trap to study recoil-induced resonances 
(RIRs) in the transient, high-gain regime. We find that two distinct mechanisms govern the atomic 
dynamics: the finite, frequency-dependent atomic response time, and momentum-space population 
redistribution. At low input probe intensities, the residual Doppler width of the atoms, combined 
with the finite atomic response time, result in a linear, transient hysteretic efTect that modifies 
the locations, widths, and magnitudes of the resulting gain spectra depending on the sign of the 
scan chirp. When larger intensities (i.e., greater than a few /xW/cm^) are incident on the atomic 
sample for several /is, hole-burning in the atomic sample's momentum distribution leads to a coherent 
population redistribution that persists for approximately 100 ^s. We propose using RIRs to engineer 
the atomic momentum distribution to enhance the nonlinear atom-photon coupfing. We present a 
numerical model, and compare the calculated and experimental results to verify our interpretation. 



I. INTRODUCTION 

Much recent work has focused on the realization of 
nonlinear optical interactions with few photons for ap- 
plication to creating more efficient all-optical components 
and to quantum information and communication schemes 
P) S 0] ■ order to drive a material into the nonlinear 
regime with only a small number of photons, the non- 
linear material must interact strongly with the incident 
radiation. An optically thick atomic sample can pro- 
vide strong atom-photon coupling, but the sample must 
be prepared in such a way that the deleterious effects 
of linear absorption are mitigated. To date, most of the 
techniques that have been proposed to increase the atom- 
photon interaction strength {e.g., EIT, cavity QED) only 
rely on manipulating internal- or cavity-atomic states. 
By combining previously-used techniques with new ap- 
proaches that use control of both internal and center- 
of-mass atomic states, novel methods for realizing en- 
hanced nonlinear optical interactions can be achieved [3] . 
Furthermore, because these methods primarily involve 
center-of-mass atomic states, they are widely-applicable 
to a broad range of atomic species and less sensitive to op- 
tical and magnetic field inhomogeneities than quantum- 
interference-based schemes. 

In this paper, we report on an approach that exploits 
the collective excitation of a spatially extended, optically 
thick sample of cold atoms. Specifically, we focus on 
using a phenomenon known as recoil-induced resonance 
(RIR), which can be viewed in terms of Raman transi- 
tions between the quantized momentum states of an atom 
[^,0]. For a given pump-probe detuning, this two-photon 
process coherently transfers atoms from one resonant mo- 
mentum state to another, thus coupling the internal and 
external center-of-mass atomic states. The atom-photon 
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coupling strength directly depends on the populations of 
atoms in the initial and final momentum states. Thus, 
by selectively engineering the instantaneous atomic mo- 
mentum distribution we can enhance the coupling 
strength and reduce the threshold for nonlinear optical 
behavior. 

In order to understand how to construct an optimal 
momentum distribution via RIRs, we must first study 
the transient dynamics of RIRs in the high-gain regime. 
For weak optical fields incident on optically thin atomic 
samples, momentum-changing RIR events alter negligi- 
bly the momentum distribution of the atoms. Thus, 
for a thermal gas of atoms with a Maxwell-Boltzmann 
momentum distribution, the resulting RIR signal has a 
Gaussian-derivative shape with small gain (loss) for neg- 
ative (positive) pump-probe detunings. For an optically 
thick atomic sample in the high-gain regime, the RIR 
feature dominates the observed spectrum Q. In order 
to understand the complex, highly- coupled dynamics in 
this case, one must consider the interplay between the 
amplification of the probe beam and modification of the 
momentum distribution as the field propagates through 
the gas of atoms. Furthermore, one must consider the 
finite response time of the material when investigating 
transient phenomena. 

As an example of self-enforced momentum engineering, 
Vengalattore et al. claim to observe enhanced nonlinear 
optical effects at low light levels giving rise to transient 
optical bistability Q. Here, we report a similar, transient 
hysteretic effect at low light levels (< 10 /iW/cm^), but 
interpret it as a purely linear phenomena. Also, we ob- 
serve a substantial modification of the atomic momentum 
distribution at higher intensities, thus demonstrating the 
feasibility of momentum-state engineering via RIRs. We 
compare the experimentally-obtained results with a nu- 
merical model to verify our interpretation of the obser- 
vations. 

This paper is organized as follows. In Section [Til we 
briefiy describe the experimental setup and Section lllll 
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describes the model used to describe the RIRs. We 
present and discuss the results in Section [IV] and Sec- 
tion |Vl respectively. Section IVll concludes the paper and 
indicates future research directions. 



II. EXPERIMENTAL SETUP 

In the present study, we create a sample of cold atoms 
via an anisotropic magneto-optic trap (MOT) as de- 
scribed in previous work ^IQ] . The sample consists of ap- 
proximately 5 X 10^ ^^Rb atoms confined in a cylindrical 
volume with a 1/e radius and length of 300 /im and 3 cm, 
respectively. A pair of counterpropagating lin_Llin laser 
beams (cooling beams) intersect the trapping volume (as 
defined by the magnetic field) at a small angle {9 ~ 10°) 
to cool the atoms along the long dimension of the trap; we 
achieve typical atomic temperatures of 20 — 30 /iK with 
this scheme. These cooling beams also act as the pump 
beams in the RIR scheme, and have an effective Rabi 
frequency of f^i/F = 2.5 (where T/2t: = 6 MHz is the 
natural linewidth of the ^Si/2{F 1) ^ ^Pz/2{F = 2) 
transition). The experiment is run in steady-state (with 
the MOT beams on), and all of the beams are typically 
detuned A = 3 — 5 F below the atomic resonance. This 
configuration enables the production of optical densities 
(OD; lout/Iin = exp(-OD)) of up to ~ 60. We con- 
trol the OD by varying the detuning of the repump laser 
beam, which allows us to investigate RIRs in both the 
low- and high-gain regimes. 

The probe beam used for RIR spectroscopy is split 
off from the pump beam, which provides the phase co- 
herence necessary for studying spcctroscopically narrow, 
multiphoton resonances. As shown in Figured] a), the 
probe beam propagates along the long axis of the trap, 
with a polarization that is parallel to that of the counter- 
propagating pump beam {i.e., the K_L configuration de- 
scribed in Ref. |ll|).The frequency detuning of the probe 
beam relative to the pump beam and spectral scan rate of 
the probe beam are independently controlled via acousto- 
optic modulators. The input probe intensity {Im) ranges 
from 0.5-200 /xW/cm^, and the scan rates vary between 
0.08-10 MHz/ms. The size of the incident probe beam 
is 100 /im (1/e^ intensity radius) with a Rayleigh range 
longer than the trap length. Although independent ex- 
periments by our group 12] and others have demon- 
strated waveguiding and focusing/defocusing effects that 
play additional roles in light propagating through such a 
dispersive atomic medium, we do not explicitly consider 
such effects here. 



III. THEORETICAL DESCRIPTION 

Most formulations of RIRs consider the atom-photon 
interaction as a stationary process, where the predicted 
spectrum is independent of the prior history of the field 
or atomic ensemble This approach is appropriate for 



© 



Probe 



© 




Eump MOT 



FIG. 1: a) Experimental RIR beam geometry. A pair of coun- 
terpropagating, lin_Llin beams, oriented at an angle 6 with 
respect to the long trap direction, act as cooling beams for 
the MOT and as pump beams in the RIR scheme. The probe 
beam propagates along the length of the trap, and is copolar- 
ized with the counterpropagating pump beam, b) RIR energy 
level scheme. The horizontal and vertical axes correspond 
to the atomic momentum and energy, respectively, and the 
parabola describes the quantized ladder of states that sat- 
isfy energy and momentum conservation. The circles repre- 
sent the relative populations of the two resonant momentum 
states. 



weak incident fields and slow probe frequency scan rates 
(relative to the atomic decay times), but one must solve 
the transient problem in the case of fast scan rates or 
substantial probe beam amplification. Theoretical [Tsj 
and experimental [l6[ studies of RIRs in the transient 
regime have been carried out in the low-gain limit for op- 
tically thin atomic samples. While the majority of studies 
focusing on high-gain center-of-mass-mediated phenom- 
ena have focused on the collective atomic recoil lasing 
(CARL) regime [13, [13 ; some recent work has also fo- 
cused on the RIR regime @. 

In order to describe the present experimental situation, 
we use a model that describes the interaction of classi- 
cal optical fields with a sample of thermal, cold atoms 
with quantized momentum states [l9| . Experimentally, 
we find that the main contribution to the RIR spectrum 
comes from the probe beam and nearly counterpropagat- 
ing pump beam (for the polarization configuration shown 
in Fig. [1]), and hence we only consider these two beams in 
the theory. Furthermore, because the angle between the 
pump and probe beams is small and the atoms are tightly 
confined in the radial direction by the trapping potential, 
we consider only motion along the longitudinal direction 
of the trap. Finally, we consider that the atoms have two 
internal states (a ground state | g) and excited state | e)) 
coupled to a quantized ladder of momentum states (see 

Fig.mb). 

We write the relevant Hamiltonian as uM 
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+ih J2 {gjay^'clik - k,)Ce{k) - H.C.)] , (1) 

where m is the atomic mass and wq is the natural fre- 
quency of the two-level atomic transition. The atom- 
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photon coupHng constant for the pump {ki,uji) and 
probe (^2 ~ —ki,uj2 = uji — S) beams are given by 
<?i,2, where gj — iJ,j[ckj / {2heoV)]^^'^ , fij is the dipole ma- 
trix element, and V is the quantization (trap) volume. 
The unitless, single-photon field amplitudes of the pump 
(probe) fields are oi (02), and c*g{k) {ce.g{k)) are cre- 
ation (annihilation) operators for the ground and excited 
states with atomic momentum hk, respectively. 

In the limit that the pump and probe beams are far- 
detuned from the atomic resonance, we can adiabatically 
eliminate the excited states. Assuming that the pump 
beam propagates with negligible attenuation, one can de- 
rive an expression for the coherence between any pair of 
initial and final momenta. By considering only the pop- 
ulations (lip = p{p,p)) of and the first-order coherences 
{rip = p{p + liP)e~*^*) between momentum states, the 
atomic evolution can be described as 



A. Thermal-Equilibrium Limit (TEL) 

We first consider the situation of arbitrary input probe 
beam powers, but we fix the atomic momentum distribu- 
tion at its thermal equilibrium value {Ilp{t) = Hth^p)- In 
this limit, excitations between the atoms (through the co- 
herence r]p) and the field (through 02) can be exchanged. 
This limit allows us to distinguish between effects pro- 
duced by population redistribution and those due to ei- 
ther the finite response time of the material or the fact 
that the atomic sample consists of an inhomogeneously- 
broadened group of radiators. Because all of these mech- 
anisms can give rise to similar transient effects, we com- 
pare both the numerical results from the full set of 
evolution equations and the fixed-population, thermal- 
equilibrium limit equations with the experimental results 
in Sec. [IVl 



tip = [-i(3*a2{-r]p + r]p-i) + c.c.]--fpop{Ilp~Ilth,p), 
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where jpop (sfcoh) are the population (coherence) decay 
rates, Ilt^.p is the thermal population distribution (typ- 
ically given by a Maxwell-Boltzmann distribution), and 
LOr ~ hkf/2'm is the single-photon recoil frequency. The 
dimensionless momentum is given hy p = ?ifc/(2?ifco) (for 
2fco = \ki — ^2!), and f3 = gig2ai/ A, where A = a;2 — wo 
is the pump-bare atomic resonance detuning (see Fig. [1] 
b). 

In order to present a self-consistent picture of the 
atom-field interaction. Maxwell's equations must also be 
solved simultaneously. Ignoring propagation effects (i.e., 
using a mean-field approximation), the time dependent 
probe field can be written as: 
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(3) 



where N is the number of atoms in the probe beam vol- 
ume, K = c/L is the free space decay of photons from 
the atomic sample of length L, ain is the amplitude of 
the input probe beam, and the summation runs over all 
momentum states. 

We numerically integrate Eqs. [211S1 using a spacing be- 
tween momentum states of 10~^p (we have verified that 
the results do not change when we use smaller steps) , and 
then sum over momenta in the range p = [—35, 35]. The 
gain experienced by the probe beam propagating through 
a distance Az is calculated as 



B. Perturbative Limit 

In order to connect our results with prior work [l5| , 
we investigate Eqs. [5][2]in the perturbative limit. Here, 
we consider only weak beams and low gain, where the 
probe amplitude and amplification of the probe beam are 
assumed to be small {i.e. 02 (i) ^ a.^ << ai). Also, as 
in the TEL, we assume that the momentum distribution 
remains at its thermal equilibrium value {Tlp{t) ~ Hth,p)- 
Under these approximations, Eqs. [2113] reduce to 



rip = i[f{p) - 5{t)]rjp - 
where f{p) = Aujr{p'^ — {p 
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1)2) and ABp = Uth,p+i - 
If 6{t) — 60 + Rt, where <5o is the probe detuning 
ai t = and R is the probe scan rate, then rip{t) can be 
explicitly solved for as 
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^M\j'^{5eff 



Rt))], (6) 



where erfi = — ierf(iz) is the imaginary error function and 
5eff = 5q — f{p) — ijcoh is the effective, complex probe 
detuning from the momentum class p. Equation [6l can 
describe both the transient behavior for a given detuning 
(for R = 0) and the RIR spectra obtained by scanning 
the probe beam (i? 7^ 0). While we do not rely on this 
model for quantitative comparison between theory and 
experiment, having an analytical solution allows us to 
gain insight into the physical mechanisms involved. 



lout/Iin = exp[2((a2 - a^n)/ain){Az/L)]. 



(4) 



This formulation allows us to model the rich, coupled 
dynamics that we describe in Sec. IIVI Before inves- 
tigating the full solution to Eqs. [UlSl we briefly discuss 
two limiting cases: the thermal-equilibrium limit and the 
perturbative limit. 



IV. RESULTS 



A. Transient Dynamics 



In order to understand the interplay between the fi- 
nite response time of the atoms and the redistribution 
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of atomic population among momentum states, we study 
RIRs over a range of input probe intensities and probe 
frequency scan rates. We begin by investigating the tem- 
poral evolution of the probe beam amplitude for a fixed 
pump-probe detuning. We note that, for the case of 
i? = in the perturbative regime (i.e., constant probe 
frequency), Eq. [H] reduces to the result for the RIR tran- 
sition rate per atom with momentum p [dPp/dt) obtained 
by Guibal et al. [H] 



dP 

oc 3[77j = (An,) 



sin((5e//t) 



(7) 



where Pp is the probability of a Raman transition be- 
tween a state with momentum p and p + 2fco- After inte- 
grating over all momentum classes, the perturbative limit 
predicts that the probe amplitude initially increases (de- 
creases) for a negative (positive) detuning, and that it 
oscillates at a frequency of S/tt before decaying to its 
steady-state value. One can interpret this result physi- 
cally by considering the individual momentum classes as 
independent oscillators. As Eq. [7| shows, each momen- 
tum class has a characteristic coherence oscillation time 
that depends on the pump-probe detuning. Thus, the 
probe field oscillations arise due to the interference of 
the radiation emitted by the various momentum classes, 
and the decay occurs as a result of net destructive inter- 
ference as the oscillators dephase relative to one another. 



b) 



;=-280 kHz 




(0 



? 0.5 



10 15 20 25 

t(|lS) 




0.5 2 5 10 

7i/5 {[is) 



FIG. 2: a) Transient RIR oscillations observed when the probe 
beam is turned on. Experimental data (solid lines) is com- 
pared with the numerical results (dashed lines) of the full 
model for S = —280 and —455 kHz. We normalize the magni- 
tude of the experimental signal to facilitate comparison with 
theory, b) Time at which the first probe maxima occurs as 
a function of the inverse detuning. Other experimental pa- 



rameters are: 
f3 = 13 kHz. 



lin = 100 ^W/cm^ T ^ 20 ^IK, N ^ 8 X lO" 



To investigate the transient regime experimentally, we 
first produce a cloud of atoms, and then rapidly turn 
on the probe beam and measure its intensity after it 
passes through the atoms (lout)- Figure [5] a) shows the 
experimentally-observed probe intensity along with nu- 
merical calculations from the full equations for two dif- 
ferent detunings. We note that, while the oscillation fre- 
quency depends only on the pump-probe detuning, the 
decay time is sensitive to the probe power, average atomic 
temperature, and the coherence decay rate. By inde- 
pendently measuring the atomic temperature (via in-situ 
RIR velocimetry [20!|) and probe power, we determine 



7co/i/27r = 20 kHz for a variety of pump-probe detunings 
by fitting the data to the full model. 

In general, we see good agreement between the numer- 
ical and experimental results over a range of detunings 
and probe powers. Figure[2]b) shows that, at short times, 
the first maxima occur at a time T^ax , which is approxi- 
mately equal to n/S, as predicted by all three models. For 
longer times, the experimentally-observed values of Tmax 
diverge from the perturbative solution. The results of the 
TEL and full-equation solutions, though, agree through- 
out all investigated times, thus indicating that, over such 
short times (several fj,s), population redistribution does 
not affect the dynamics. 
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FIG. 3: a) Experimentally-observed and numerically- 
predicted temporal evolution of the probe intensity occuring 
after the probe intensity is switched from 200 to 1 /iW/cm^. 
b) Relative hole depth resulting from different initial probe 
beam intensities. Other experimental conditions are: T = 20 
/iK, iV = 8 X 10^ (3 = 13 kHz, and S = -140 kHz. 

At slightly longer times, the effects of population re- 
distribution and rethermalization become important. In 
order to study this regime, we performed an experiment 
similar to Ref. @. We illuminate the atoms with an in- 
tense probe beam for ~ 100 ms before rapidly dropping 
the probe beam intensity to 1 /zW/cm^. For large ini- 
tial intensities, we measure small, weak-beam gains that 
slowly increase to their steady-state values. To prop- 
erly interpret this response, we need to distinguish be- 
tween two different effects: the initial, rapid oscillations 
in the output probe intensity (as described in the previ- 
ous paragraph) and the slower response due to changes in 
the momentum-space populations. Figure |3| shows that, 
even when the momentum distribution is fixed, we ex- 
pect a small decrease in gain after we reduce the probe 
intensity. This occurs because of the destructive inter- 
ference between adjacent momentum classes excited by 
the strong beam that relax at different rates. In agree- 
ment with Vengalattore et ai, we interpret the slower 
relaxation to steady state as due to a redistribution of 
atomic population. The intense beam removes atomic 
population from a particular momentum class, thereby 
reducing the gain; once the probe intensity is reduced, 
random scattering and collisions rethermalize the atoms, 
thus refilling the hole burned in the momentum distribu- 
tion and increasing the gain. 

We estimate the time necessary to remove a substan- 
tial amount of population for given an incident power as 
tpump ~ 27r/rii2 (where fli2 = 4/3a2 is the two-photon 
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Rabi frequency). For typical experimental parameters 
{i.e., P — IS kHz and a beam diameter of 100 /im), the 
time required to redistribute population is given as 

tpump[^^s] - 167 X (/„[/iW/cni'])-i/2_ (8) 

Equation [8] predicts a population redistribution time 
of tpump — 53 fis for an intensity of /i„ = 10 
/xW/cm^, corresponding to an input energy density of 
~2 photons/(A^/27r). In the following paragraphs, we 
compare these predictions with the numerical and exper- 
imental results. 

Figure [3] a) shows a representative plot of the normal- 
ized probe gain {lout/Iss, where Igs is the steady-state 
probe intensity) as a function of time for an initial probe 
intensity of 200 /iW/cm^. By fitting the measured signal 
with the full experimental model, we extract a rether- 
malization rate of 7pop/27r = 3.4 kHz. We repeat this 
experiment for a range of initial probe intensities, and 
observe that the relative hole depth {dhole = Imin/Iss, 
where Imin is the minimum transient gain) decreases with 
the initial intensity. Furthermore, dhoie is equal to unity 
for intensities less than ~ 10 /zW/cm^, indicating that 
momentum-space population redistribution negligibly ef- 
fects the atom-field dynamics for sufficiently low intensi- 
ties. 

To verify our interpretation of the temporal evolution 
of the gain, we compare the numerical results from the 
full and TEL calculations with the experimental data. 
Figure [3] b) shows the values of dhole as a function of the 
strong probe beam intensity. The TEL predicts a slight 
reduction in gain due to transient effects identical to 
those described above and does not agree with the data. 
On the other hand, the full solutions accurately predict 
both the relative hole depth and decay time. Thus, in 
agreement with Eq. [51 we find that population redistri- 
bution plays an important role when intensities greater 
than a few /iW/cm^ are incident upon the atoms for tens 
of /is. 

B. Transient Hysteresis 

We also investigate the effects of momentum redistri- 
bution by studying the RIR spectrum. We measure the 
spectrum by scanning the probe beam frequency across 
the resonance at a fixed rate R. While, in steady state, 
the resulting RIR signal reflects the derivative of the 
equilibrium momentum distribution, scan rates that tra- 
verse the resonance on time scales comparable to the RIR 
turn-on dynamics (~ 1 fis) produce history-dependent 
spectra that reflect the local momentum distribution. In 
this section, we focus on the effects of the frequency- 
dependent, finite response time of the material and the 
controllable redistribution of momentum-space popula- 
tion on the RIR spectra. 

Figure [?] shows the RIR spectra for a low input probe 
intensity (0.5 /iW/cm^) at different scan rates. For slow 
scans, positive {R > 0) and negative {R < 0) chirps 
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FIG. 4: a), c) Experimental and b), d) theoretical (T — 20 
fiK, iV = 8 X 10'', /3 = 13 kHz) RIR spectra at slow and fast 
scan rates, respectively, for a probe intensity of 0.5 /xW/cm^. 
Solid (dashed) lines correspond to negative (positive) chirps. 



produce identical spectra, whereas fast scans result in 
resonance line shapes with altered locations of the max- 
imum value of the gain, widths, and amplitudes. Figure 
|4] a) and c) show experimental spectra obtained with a 
slow (i? = 0.1 MHz/ms) and fast (i? = 8 MHz/ms) scan 
rate, respectively. The slow scan probes the steady-state 
response of the sample, and we observe no chirp depen- 
dence. For a fast scan with a negative chirp, though, the 
spectrum has a larger gain peak, narrower gain feature, 
and a larger shift from 6 = relative to the steady-state 
case. For a positive chirp, the opposite is true. The chirp- 
dependence of the numerically-calculated spectra, shown 
in Fig. [3] b) and d) , agree well with the experimental 
data. We note, though, that the amplitude of the gains 
predicted by the model do not match the experimentally- 
measured values. We believe that this discrepancy is 
caused by two main effects not accounted for in the 
model: additional resonances caused by the MOT beams 
(producing small, additional Raman resonances between 
light-shifted levels [2l|, [S^l), and pr opagati on effects that 
the mean-field assumption ignores [ij, [l^ . 

One can understand this transient hysteresis in terms 
of the short-time dynamics of the atomic sample. We 
first note that, because the momentum distribution is 
not significantly modified at very low probe powers ap- 
plied for short times (see Sec. IIV A[) . one cannot explain 
the results of Fig. [4] in terms of population redistri- 
bution. Rather, the hysteretic effect occurs as a result 
of the frequency-dependent response time of the atoms. 
As mentioned above, the time it takes for a particular 
momentum class to reach its first maximum is approxi- 
mately inversely proportional to the pump-probe detun- 
ing. When the probe frequency is scanned farther from 
(5 = 0, the resonant momentum classes respond increas- 
ingly quickly as 6 increases, resulting in a situation where 
a range of momentum classes radiate at their maximum 
intensity at the same time. In the opposite case, when the 
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probe frequency is scanned toward 6 = 0, the resonant 
momentum classes reach their maximum radiated pow- 
ers at different times. Considering the resulting probe 
intensity as a superposition of this inhomogeneous collec- 
tion of radiators thus explains the increased (decreased) 
and narrowed (broadened) gain peak for the fast, neg- 
atively (positively) chirped case. In a similar way, the 
finite response time of the material effectively delays the 
occurrence of the gain peak for either case of the chirp, 
resulting in the observed shifts in the gain peaks. 
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FIG. 5: a) and c) Experimental results for the RIR peak gain 
as a function of scan rate for an input intensity of 0.5 and 
100 fiW/cia^, respectively. Boxes and circles correspond to 
positive and negative chirps, b) and d) Theoretical results 
for the RIR peak gain as a function of scan rate for an input 
intensity of 0.5 and 100 /xW/cm^, respectively. Solid (dashed) 
lines correspond to solutions to the full (TEL) equations for 
positive and negative chirps, as indicated. Other experimen- 
tal conditions are: T = 20 ^K, iV = 8 x 10** and 13 = 13 
kHz. 

We confirm our interpretation of the hysteresis at low 
powers by first studying the perturbative solution given 
in Eq. [6l For finite scan rates, Eq. [6] produces the 
experimentally-observed, chirp-dependent variations of 
the RIR gain feature. As this model explicitly ignores 
changes to the momentum distribution as well as any 
back-action between the atoms and photons, one can un- 
derstand the observed hysteresis as a transient, linear ef- 
fect. This result is directly analogous to the modification 
of the resonance lineshape of a damped, driven harmonic 
oscillator when the driving frequency is scanned across 
the resonance. 

To quantify the importance of population redistribu- 
tion in the experimentally-observed spectra, we compare 
the results of the full and TEL calculations. Figure [H] 
b) shows that both calculations predict almost the same 
chirp-dependence of the gain, thus indicating that mo- 
mentum redistribution does not significantly contribute 
to the observed gain. Figure [5] a) shows experimental 
data that demonstrates that, for slow scans (i? < 0.5 
MHz/ms), the atoms reach steady-state and the chirp- 
dependence disappears. At faster scan rates, the gain 
increases, as discussed above. Beyond i? ~ 3 MHz/ms, 



though, the gain observed for both chirps decreases as 
the atom-photon interaction time becomes too short. 

At higher powers, population redistribution plays a 
role in the observed RIR spectra. Figure[5]c) and d) show 
the measured and calculated dependence of the peak gain 
value on the scan rate for an input probe intensity of 100 
/iW/cm^. For fast scan rates, the numerical results of 
the thermal and full model coincide, thus indicating that 
the probe beam does not spend long enough at each de- 
tuning to significantly modify the populations. For slow 
scan rates, the results of the two models diverge. In this 
region, momentum redistribution decreases the predicted 
gain by effectively reducing the population difference be- 
tween the resonant momentum states. This situation dif- 
fers from Ref. [§] because the scan rate considered here 
is slow enough to allow for new, equilibrium momentum 
distributions to occur for each resonant momentum class 




0.1 0.51 



510 



0.1 0.51 



5 10 
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FIG. 6: a) RIR peak gain and b) shift of the peak gain loca- 
tion as a function of scan rate. Boxes and circles correspond 
to experimental results, and the solid and dashed lines corre- 
spond to theoretical results for input intensities of 0.5 and 100 
/iW/cm^, respectively. Other experimental conditions are: 
T = 20 mK, = 8 X 10^ [3 = 13 kHz, and S = -140 kHz. 



The relative differences between the peak gain magni- 
tude and location help further clarify the situation. Fig- 
ure [S] a) and b) show that the ratio of the negative (g-) 
and positive (5+) chirp peak gains and the relative shift 
of the location of the gain peak {Av^ax = i5™°^ — (5™°^) 
agree qualitatively with results of the full model and do 
not strongly depend on intensity. Thus, while popu- 
lation redistribution between momentum classes occurs 
in this system, the time scales over which it acts are 
longer than the time the probe spends resonant with 
a given momentum class for fast scans. We therefore 
conclude that the observed transient hysteresis results 
from the linear, frequency-dependent response time of 
the inhomogeneously-broadened system, rather than co- 
herent population redistribution. 



V. DISCUSSION 

Our results are consistent with those described by Ven- 
galattore et al. Q; only our interpretation of the ob- 
served phenomenon differs. For the mechanism described 
by Vengalattore et al., a negatively-chirped scan shuffles 
population such that the population in the ground mo- 
mentum state increases, thus increasing the population 
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difference between the resonant momentum states and 
enhancing the gain. For a positive chirp, though, one 
should observe a similar enhancement as the population 
difference is similarly increased by a removal of popula- 
tion from the excited momentum state. As we observe 
only an increase in gain for a negative chirp, in agree- 
ment with the predictions of the perturbative model, we 
conclude that the effect described in Ref. @ does not 
play a major role in our system. Nevertheless, we clearly 
demonstrate effects of momentum redistribution. 



VI. CONCLUSION 

In conclusion, we have investigated RIRs in the tran- 
sient, high-gain regime, which supplements the work de- 
scribed in Ref. By studying the RIR signal pro- 
duced at short times for fixed frequencies, we measure 
the effective population and coherence decay rates. Also, 
we note two important effects that influence the RIR 
spectrum: the frequency-dependent, finite response time 
of the material, and the redistribution of momentum- 
space population for sufficiently high probe intensities. 



By measuring the RIR spectrum for various powers and 
probe frequency scan rates, we observe transient hys- 
teretic phenomena that arise from both linear (finite ma- 
terial response time) and nonlinear (population redistri- 
bution) effects. The results of this study demonstrate 
that momentum-space hole-burning is possible because 
of the inhomogeneously-broadened nature of the RIR, 
and that it persists for approximately 100 /is. By tai- 
loring the atomic momentum distribution via coherent 
population redistribution mediated by RIRs, the nonlin- 
ear atom-photon coupling can be controlled to enhance 
or reduce the nonlinearity. Also, optimizing the setup 
to make use of collective effects resulting from atomic 
bunching in position space can further improve the ob- 
served nonlinear coupling. Together, these effects make 
this system an excellent candidate for the realization of 
ultra- low- light nonlinear optics. 
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